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Abstract 

We report measurements of the charm-anticharm production asymmetries for A+, 
TP C , Sc*, and A+(2625) baryons from the Fermilab photoproduction 

experiment FOCUS (E831). These asymmetries are integrated over the region where 
the spectrometer has good acceptance. In addition, we have obtained results for 
the photoproduction asymmetries of the A c baryons as functions of pl, p\, and 
xf- The integrated asymmetry for A+ production, (c A + — cr A -) / (o" A + + c A -), 
is 0.111 ± 0.018 ± 0.012, significantly different from zero. The asymmetries of the 
excited states are consistent with the A c asymmetry. 



The FOCUS experiment uses a photon beam on a Beryllium Oxide target to 
produce charm particles. In high energy photon-hadron interactions, pairs of 
charm-anticharm quarks are produced predominantly through photon-gluon 
fusion [1]. At leading order in Quantum Chromodynamics (QCD), the pro- 
duced charm and anticharm particles are identically distributed in the kine- 
matic variables. At next-to-leading order, small asymmetries are expected 
between charm and anticharm production [2,3,4]. However, these predicted 
asymmetries would be too small to be measured at the current level of exper- 
imental statistics. Interactions between the produced charm quarks and those 
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in the struck nucleoli during hadronization can also induce production asym- 
metries. One common model is that of string fragmentation as implemented 
in PYTHIA [5]. In this model, the charm and anticharm quarks are connected 
through color strings to the quarks in the struck nucleon. The energy in these 
strings is converted to particles by "popping" qq pairs out of the vacuum. The 
simplest asymmetry example occurs when a charm quark is connected to the 
diquark in the nucleon by a low energy string with insufficient energy to pro- 
duce any additional particles. In this case, the charm quark can combine with 
valence u and d quarks to form a A c while a c quark can only form mesons 
when it combines with the valence quarks. 

In this letter we present a high-statistics measurement of the production asym- 
metry of A c baryons from photon-nucleon interactions, providing the first con- 
vincing evidence for a non-zero asymmetry. This letter also contains the first 
measurements of this asymmetry as functions of pl, v\i an d %f- In addition, 
the production asymmetries of the excited charm baryons E°*, 
and A+(2625) are presented for the first time. All of these states decay to a 
A c plus one or two charged pions. 

The FOCUS (Fermilab E831) experiment was designed to study charm parti- 
cle physics. Charmed hadrons are produced by the interaction of high energy 
photons ((E) ~ 175 GeV for events in which a charm decay was reconstructed) 
with a segmented Beryllium Oxide (BeO) target. The photons are produced 
by bremsstrahlung in a lead target with a 300 GeV e + /e~ beam. Vertex recon- 
struction is performed using four silicon strip planes interleaved with segments 
of the target followed by a 12 plane silicon strip vertex detector. Downstream 
of the vertex detector, tracking and momentum measurements are made using 
five stations of multiwire proportional chambers and two large aperture mag- 
nets with opposite polarity. Three multicell Cerenkov counters operating in 
threshold mode are used to identify electrons, pions, kaons, and protons over 
a wide range of momenta. The spectrometer also contains a hadron calorime- 
ter, two electromagnetic calorimeters, and two muon detectors. Data were 
collected during the 1996-97 fixed-target run. 

The A+ particles are reconstructed using the pK~n + decay mode. 1 The decay 
vertex is formed from three charged tracks in the event. The vector sum of their 
momenta is projected back toward the target and used as a seed to intersect 
with at least one other track in the event to form a production vertex. We 
require the production and decay vertices to be separated by at least 5.5 
<j£, where <j£ is the uncertainty in the measured vertex separation. We also 
place goodness-of-fit criteria on the primary and secondary vertices, requiring 
that the confidence level for each vertex be greater than 1%. Identification of 
the decay tracks by particle type is performed using the Cerenkov detector 
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system [6]. A x 2 -like variable is formed using the on/off status of all cells 
within a particle's Cerenkov cone {(3=1). For each of the four possibilities, 
electron, pion, kaon, and proton, we calculate Wi = —2 S^ ens log Pj, where i 

is the particle type and j the cell number. Pj is the probability that the j 
cell will yield the observed response given particle type i. Identification is then 
based on differences in the W^. The proton candidate is required to have the 
proton hypothesis favored over the kaon and pion hypotheses by 1 and 4 units 
of log likelihood, respectively. The kaon hypothesis for the kaon candidate is 
required to be favored over the pion hypothesis by 3 such units. For the pion 
candidate, the pion hypothesis cannot be disfavored by more than 6 units of log 
likelihood relative to the most likely hypothesis. This very loose requirement is 
also applied to the pions from the excited charm decays, described below. To 
remove longer lived charm backgrounds, the lifetime of the A c in its rest frame 
must be shorter than 8 times the world average lifetime [7] . A c candidates are 
identified as all those events which satisfy the above criteria and which fall 
into the mass range between 2.10 and 2.45 GeV/c 2 . Finally, we restrict our 
sample to be in the (large) phase space region given by 40 < p L < 200 GeV/c 
and pi < 6.0 (GeV/c) 2 . 

The and candidates are reconstructed by combining the 

A+ candidates within 2a of the mean A+ mass with a single charged pion track. 
All possible combinations in the event are tried. The confidence level of the 
E c decay vertex is required to be greater than 1%. To reduce systematic errors 
coming from the reconstruction of the A+ and obtain better signal-to-noise, 
we study the S c states using the difference in the invariant mass of each E c 
state and the invariant mass of the A+, AM. Because the pion is typically of 
low momentum, it suffers from a considerable amount of multiple scattering, 
and the uncertainty on its momentum dominates the error on the S c invariant 
mass. To improve the momentum measurement, the primary vertex is refit 
without this soft pion track (if possible), and the pion direction is recomputed, 
forcing it to come from the refit primary. 

A+(2625) candidates are reconstructed by combining the A+ candidates within 
2a of the mean A+ mass with all combinations of two pions of opposite charge 
in the event. As for the S c states, we use the mass difference plots and force 
the two pions to come from the primary. 

FOCUS data were taken with three different beam energies and with two 
different radiators. Since the asymmetries can depend on the energy spectrum 
of the incident beam photons, we include only those events in our sample that 
come from data taken with a 300 GeV e~/e + beam on a lead radiator of 20% 
of a radiation length. This selection results in a charm baryon yield which is 
about 75% of the yield found using all of the data. For reconstructed charm 
events passing the trigger, the beam energy is reasonably well described by a 
Gaussian with mean of 175 GeV and width of 45 GeV. 
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From studies of high statistics meson decays we find no evidence that there is 
any asymmetry introduced by charge bias in the spectrometer. Even so, since 
the acceptance depends on the longitudinal and transverse momenta of the 
produced baryon state, if these spectra are different for particle and antipar- 
ticle over the range for which the asymmetry is measured, there will be an 
acceptance difference for the two samples which can be corrected. This cor- 
rected asymmetry gives the asymmetry for an experiment with flat acceptance 
in p L and p\ over the range 40 < p L < 200 GeV/c and p\ < 6.0 (GeV/c) 2 . 

To compensate for this acceptance difference as a function of p\ and p L , the 
production asymmetry A is determined by: 

A = (1) 
N/e + N/e V ; 



where N and iV are the numbers of reconstructed baryons and antibaryons, 
respectively and e(e) is the baryon (antibaryon) reconstruction efficiency. The 
efficiencies e and e are calculated using the FOCUS Monte Carlo simulation 
program. The detector simulation uses a detailed description of the FOCUS 
detector. The physics processes are generated using Pythia (version 6.127), 
with many Pythia parameters tuned to match the observed FOCUS physics 
distributions. The remaining discrepancy between data and Monte Carlo is 
removed by weighting Monte Carlo events to exactly match the observed data 
Pl and p\ distributions. For the excited charm states, the particle/antiparticle 
efficiencies are used to correct the asymmetry. For the higher statistics A c 
decays a further step is taken. Since very little p\ dependence on efficiency is 
observed, the Monte Carlo events are used to obtain the efficiency variation 
versus p^. This efficiency is fit to a function which is then used to weight each 
event. 

The global A c asymmetry is obtained by fitting each of the two weighted 
invariant mass data distributions shown in Fig. 1 with a Gaussian signal and 
quadratic background. The Gaussian mean and width are allowed to float 
separately for A+ and A~. The unweighted yields for A+ and A~ are 5427±120 
and 4242±108, respectively. 

To account for natural widths and changing resolutions, the excited charm 
baryon states are fit slightly differently. The Monte Carlo is used to generate 
the signal shape which is fit to a spline function. This spline function, properly 
normalized, is used as the signal shape for the data, with the mass allowed to 
float. The advantage in the case of the S c , and S* is that the spline function is 
able to account for the detector resolution and natural width of the state. For 
the A* state, the spline function is a better model of the detector resolution 
than a Gaussian due to the small phase space. The background shapes for S c , 
S* and A c (2625) are a threshold function JV(1 + a(AM -m^AM 13 ), quadratic 
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Fig. 1. pK~ir + weighted invariant mass distribution split by charge. The weights are 
used to account for efficiency loss versus momentum on an event-by-event basis. The 
average weight for the whole histogram is given in the figure as w. The unweighted 
yields are also shown. 

polynomial, and linear polynomial, respectively. The fits to each data sample 
are shown on the corresponding data plots in Figs 2-6. The global asymmetries 
are calculated from the returned yields. 



The A c data samples are divided into bins of each of the kinematic variables 
Pl, Pt, and xf, integrating over the full range of the other variables. The mass 
plots for each bin are fit using a Gaussian signal and quadratic background, 
as for the full data set. For these fits, the A+ (A~) mass is fixed to the mass 
obtained by fitting the full A+ (A~) sample. The widths are fixed to follow the 
Monte Carlo widths for each bin. The yields from these fits are used to obtain 
the production asymmetry versus pl, Pt, an d Xf- The efficiency corrected 
asymmetry distributions for A c vs pi, Pt, an d Xf are shown in Figs. 7, 8, 
and 9. The xp measurement requires knowledge of the incoming beam energy 
which is only available in about 30% of the data sample and is therefore of 
lower statistics. 

In our first systematic error check, background studies were performed to as- 
sure that feedthrough from other charm states does not contribute to the 
asymmetry. The systematic errors are obtained by making the same measure- 
ments under different analysis conditions. We performed these studies for the 



6 




Fig. 3. A+7r invariant mass distribution split by charge. 
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Fig. 6. A+7r + 7r~ invariant mass distribution split by charge, 
global asymmetry and the asymmetry versus the kinematic variables reported. 

For the A c and the three excited states (E c , E*, and A*), fits with different bin 
widths and with bins shifted by 1/2 bin were performed. A variation which 
fixed the mass for baryons and antibaryons to the value from the total sample 
was also studied. We also investigated differences due to the fit functions 
used. For the E c states a Gaussian function with a quadratic background was 
utilized to fit the signals. For the E* mass distributions we used a linear and 
cubic fit for the background instead of the quadratic background. For E c and 
E*, fits with spline functions from two additional Monte Carlo samples where 
the natural widths were varied by ±lcr were used to estimate the uncertainty 
in the knowledge of the natural widths. For the A* we used two different fit 
functions to provide a systematic check: a Gaussian for the signal shape and 
a linear background and the A* spline function with a quadratic background. 
Additionally, for all the states we used the respective raw asymmetries to 
check for systematic problems with the efficiency correction. For the A c we 
also calculated the asymmetry using the efficiencies directly from the Monte 
Carlo instead of using the Monte Carlo to obtain an efficiency function. 

The r.m.s. of all of these variations is used as an estimate of the systematic 
error. 

The global asymmetries for all of the charm baryons studied in this analysis 
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P L (A C ) (GeV/c) 

Fig. 7. pl asymmetry distribution for A c along with statistical and systematic errors 
and the Pythia prediction - the systematic errors are superposed on the statistical 
errors and are smaller than the statistical errors in every bin. 




Fig. 8. -p\ asymmetry distribution for A c along with statistical and systematic errors 
and the Pythia prediction - the systematic errors are superposed on the statistical 
errors and are smaller than the statistical errors in every bin. 

are shown in Table 1. The first errors are statistical and the second are system- 
atic. Table 1 also shows comparisons to Pythia asymmetries calculated with 
unreconstructed events. The Pythia predictions come from running version 
6.203 with all parameters left at the default setting. The events are generated 
with the correct beam energy distribution and only candidates within the 
nominal phase space, 40 < pl < 200 GeV/c and p\ < 6.0 (GeV/c) 2 , are used 
in determining the asymmetry. The E691 and E687 photoproduction experi- 
ments have previously reported global asymmetries for A c of 0.110 ± 0.089 [8] 
and 0.035 ± 0.076 [9], respectively. Our results are similar to those obtained 
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Fig. 9. xf asymmetry distribution for A c along with statistical and systematic errors 
and the Pythia prediction - the systematic errors are superposed on the statistical 
errors and are smaller than the statistical errors in every bin. 

by these experiments although comparisons are not straightforward since all 
three experiments have different phase space and beam energy distributions. 
Table 1 also shows the efficiency ratios of particles to antiparticles for the 
charm baryons studied. 

Table 1 

Raw and corrected global production asymmetry for the charm baryons compared 
to the predictions of default Pythia version 6.203. The efficiency ratio of particles 



to antiparticles is also shown. 


Baryon 


Raw Asymmetry 


Corrected Asymmetry 


Pythia 


Efficiency Ratio 


A c + 


0.123 ±0.017 


0.111 ±0.018 ±0.012 


0.073 


1.023 ± 0.005 


£++ 


0.147 ±0.073 


0.136 ±0.073 ±0.036 


0.126 


1.024 ±0.010 


s° 


0.005 ± 0.089 


0.005 ± 0.089 ± 0.024 


0.128 


1.000 ± 0.009 


£++* 


0.188 ±0.105 


0.181 ±0.105 ±0.033 


0.133 


1.016 ±0.006 


£°* 


0.299 ±0.165 


0.298 ±0.165 ±0.023 


0.132 


1.002 ± 0.006 


A+(2625) 


-0.066 ± 0.086 


-0.075 ±0.087 ±0.021 


N/A 


1.019 ±0.017 



We have studied the photoproduction asymmetry of A+ versus A~ using the 
decay channel pK~n + . From ~ 10,000 A c events we present the first results 
of this asymmetry as functions of pl, p?, and xp. These results show a clear 
positive asymmetry. The global A+ asymmetry is measured to be 0.111 ± 
0.018 ±0.012. 

The production asymmetry of excited charm states which decay to A+, in- 
cluding the £++, £°, £++*, £°*, and A+(2625) was also measured for the first 
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time. The measurements generally indicate a positive asymmetry similar to 
the A c . Because of the smaller sample size, however, they are also consistent 
with zero. 

We find that the string fragmentation model as implemented in Pythia ver- 
sion 6.203 does not describe the A+ asymmetry dependence on p^, pj,, or xp. 
Our measurements indicate that the asymmetry shows no significant depen- 
dence in these variables. 
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